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Abstract 

In this article, we present emerging fields of quantum chemistry at 
finite temperature. We discuss its recent developments on both exper- 
imental and theoretical fronts. First, we describe several experimen- 
tal investigations related to the temperature effects on the structures, 
electronic spectra, or bond rupture forces for molecules. These include 
the analysis of the temperature impact on the pathway shifts for the 
protein unfolding by atomic force microscopy (AFM), the tempera- 
ture dependence of the absorption spectra of electrons in solvents, and 
the temperature influence over the intermolecular forces measured by 
the AFM. On the theoretical side, we review advancements made by 
the author in the coming fields of quantum chemistry at finite temper- 
ature. Starting from the Block equation, we have derived the sets of 
hierarchy equations for the reduced density operators in both canon- 
ical and grand canonical ensembles. They provide a law according 
to which the reduced density operators vary in temperature for the 
identical and interacting many-body systems. By taking the indepen- 
dent particle approximation, we have solved the equations in the case 
of a grand canonical ensemble, and obtained an energy eigenequation 
for the molecular orbitals at finite temperature. The explicit expres- 
sion for the temperature-dependent Fock operator is also given. They 
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form a mathematical foundation for the examination of the molecu- 
lar electronic structures and their interplay with finite temperature. 
Moreover, we clarify the physics concerning the temperature effects 
on the electronic structures or processes of the molecules, which is 
crucial for both theoretical understanding and computation. Finally, 
we summarize our discussion and point out some of the theoretical 
and computational issues in the fields of quantum chemistry at finite 
temperature for the future exploration. 

Keywords Quantum chemistry at finite temperature; temperature depen- 
dent; polymers; protein folding; protein unfolding; intermolecular forces; sol- 
vated electrons; molecular crystals 

1 Introduction 

The history for quantum chemistry development is almost synchronous to 
that of quantum mechanics itself. It begins with Heitler and London's inves- 
tigation of the electronic structure for H 2 molecule shortly after the establish- 
ment of the wave mechanics for quantum particles pQ. There are two major 
types of molecular electronic theories: valence bond approach vs. molecular 
orbital method with the latter being the popular one for the present study. 
It has gone through the stages from the evaluation of molecular integrals 
via a semiempirical way to the one by an ab initio method. The correlation 
issue is always a bottleneck for the computational quantum chemistry and 
is under intensive examination for over fifty years [2]. For large molecular 
systems such as biomolecules and molecular materials, the development of 
the combined QM / MM approach, the pseudopotential method and the lin- 
ear scaling algorithm has significantly advanced our understanding of their 
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structure and dynamics. There are about eight Nobel prize laureates whose 
researches are related to the molecular electronic structure theory. This not 
only recognizes the most eminent scientists who have made the outstanding 
contributions to the fields of quantum chemistry, but more importantly, it 
indicates the essential roles the electronic structure theory has been play- 
ing in the theoretical chemistry as well as for the whole areas of molecular 
sciences 0. Nowadays, quantum chemistry has been becoming a maturing 
science [3J 0] . 

Nevertheless, the modern fields of quantum chemistry are only part of the 
stories for the molecular electronic structure theory. From the pedagogical 
points of view, the quantum mechanics based on which the traditional quan- 
tum chemistry is built is a special case of more general quantum statistical 
mechanics jHUHlIZI- in reality, the experimental observations are made under 
the conditions with thermodynamic constraints. Henceforth, there is a need 
to extend the contemporary areas of quantum chemistry to the realm of, for 
instance, finite temperature 0EIIZI- 

Indeed, many experimental investigations in various fields and for dif- 
ferent systems have already shown the temperature or pressure effects on 
their microscopic structures [8-31,50-58,64,72-88]. The polymeric molecule 
is one of the most interesting systems for this sort of studies [8-17]. The 
experimental measurement on the absorption spectra, photoluminescence 
{PL), and photoluminescence excitation (PLE), and spectral line narrow- 
ing (SLN) for the PPV and its derivatives all show the same trend of the 
blue shift with an increasing temperature [HI EH] . This attributes to the 
temperature dependence of their very rich intrinsic structures such as the 
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vibronic coupling jTHJEJlE!- The experimental inquiry of the temperature 
impacts on the biomolecules started in the late nineteenth century [THl HH] . 
Most recently, it has been extended to the study of folding and unfolding of 
proteins or DNAs |2TJ] El 122] • in addition to the observed patterns for the 
unfolding forces with respect to the extension or temperature, it has been 
proved that the temperature-induced unfolding is another way for the exam- 
ination of mechanisms or pathways of protein folding or unfolding processes 
[20-24]. The newest related development is on the AFM measurement made 
by Lo et al. of the intermolecular forces for the biotin-avidin system in the 
temperature range from 286 to 310K [64] . It has shown that an increase of 
temperature will almost linearly decrease the strength of the bond rupture 
force for the individual biotin-avidin pair. The analysis of the temperature 
influences over the absorption spectra of the solvated electron began in the 
1950's and it is still of current interest. A striking effect is that an increasing 
temperature will cause the positions of their maximal absorption red shift 
[72-86]. 

In two papers recently published, we have deduced an energy eigenequa- 
tion for the molecular orbitals [3 EJ. It is the extension from the usual 
Hartree — Fock equation at zero temperature to the one at any finite temper- 
ature [SHI ED] • It opens an avenue for the study of the temperature impacts 
on the electronic structures as well as their interplay with the thermody- 
namic properties. In the third section, we will present this equation and give 
the details for its derivation. We will also expound the physics concerning 
the effects from temperature and classify them into two types. The one is 
at the single orbital level and the other is really an iV-body effect. In the 
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next section, we will show four major types of experiments related to the 
investigation of the temperature influences over the microscopic structures 
of molecular systems. In the final section, we will discuss and analyze our 
presentations, and point out both theoretical and computational issues for 
the future examination. 



2 Experimental Development 

In this section, we focus our description of the experimental investigations 
related to the temperature effects on the bondings, structures, or electronic 
spectra of molecules. We choose four kinds of the most recent developments 
in these fields which are of chemical, biological, or material interests. 

2.1 Temperature effects on geometric structure and 
UV-visible electronic spectra of polymers 

The first important systems where the major issues related to the temper- 
ature influences over the geometric structures and electronic spectra are 
the polymeric molecules. Many experimental inquiries and some theoreti- 
cal works already exist in the literature [8-17]. However, how temperature 
changes the microscopic structures of the polymers are still not completely 
understood and there are many unresolved problems in interpreting their 
electronic spectra. We list here a few very interesting experimental exami- 
nations for the purpose of demonstration. 

The poly(^-phenylenevinylene)(PPV A ) is one of the prototype polymeric 
systems for the study of their various mechanical, electronic, and optical 
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properties. The impact from temperature on the absorption spectra, the 
photoluminescence (PL), and the photoluminescence excitation (PLE) of 
the PPV has also been investigated both experimentally and theoretically [HI 
In an experiment performed by Yu et al., the absorption spectra are 
measured for the PPV sample from the temperature 10 to 330-fT. The details 
of the experiment are given in their paper [10J. The resulting spectra for the 
absorption at T = 80 and 300-fT are shown in Figure 1 of that paper. We 
see that there is a pronounced change in the spectra when increasing the 
temperature. They also study the PL and PLE spectra for the PPV . The 
measured PL spectra at two temperatures: 77 and 3007^ are demonstrated 
in Figure 3, and the PLE spectra at those temperatures are depicted in 
Figure 4 of the paper [lUj . They both show the dramatic changes of the 
band blue shift when the temperature is increased. Similar studies have also 
been carried out before by the other groups [HI E] • They observed the similar 
behaviors. 

Another interesting analysis is related to the temperature effects on the 
spectral line narrowing (SLN) of the poly(2-methoxy-5-(2'-ethylhexyloxy)- 
1,4-phenylenevinylene) (MEH — PPV) spin-coated from either THF or CB 
solvents ^2]- I n the experiment conducted by Sheridan et al., the SLN is 
measured together with the absorption and PL as shown in Fig. 1 of their 
paper. It is found that the same trend of the SLN blue shift is observed 
as that for the absorption and PL with an increasing temperature. They 
attribute this to the same reason of the electronic structure modification 
resulting from the variation of temperature. 



6 



2.2 Temperature effects on structure, dynamics, and 
folding/unfolding of biomolecules 

Biomolecules are complex systems, featuring a large molecular size, a hetero- 
geneity of atomic constitutes and a variety of conformations or configurations. 
Their energy landscapes thereby exhibit multiple substates and multiple en- 
ergy barriers, and vary in size for the barrier heights [213 1213 123 12H1 l2Hj • The 
temperature should have a strong influence over their structure and dynamics 
including the folding or unfolding [18-58]. This effect could be either from the 
fluctuation of thermal motions of the molecules or due to the redistribution 
of the electronic charge as we will discuss in the next section. 

The experimental observation of the temperature impact on the micro- 
scopic structure of biological systems dates back to the very early days. One 
focus, for example, is on the measurement of the elastic properties of the hu- 
man red blood cell membrane as a function of temperature ^HlCSl- Another 
related investigation is about the influence over the thermal structural tran- 
sition of the young or unfractionated red blood cells due to the involvement 
of the protein spectrin which might modify the spectrin-membrane interac- 
tion [SHI EH]. Most recently, the atomic force microscopy (AFM) has been 
used to detect the impact from the variation of temperature on the spectrin 
protein unfolding force as well as on the bond rupture force for the biotin- 
avidin system |2ffl I2H 1221 EH ■ 

The AFM is a surface imaging technique with an atomic-scale resolu- 
tion capable of measuring any types of the forces as small as 10 -18 N. It 
combines the principles of the scanning tunnelling microscopes (STM) and 
the stylus profilometer, and therefore can probe the surfaces of both con- 
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ducting and nonconducting samples [321 EH] • The imaging on soft materials 
such as biomolecules with the AFM has been accomplished beginning in the 
1980's jSH ESI EH]- Recently, it has been applied to measure the adhesive 
forces and energies between the biotin and avidin pair as we will show in the 
next subsection [37J EHl EH HO] ■ Unlike other experimental techniques, the 
AFM features a high precision and sensitivity to probe the surface with a 
molecular resolution, and can be done in physiological environments. 

In an AFM investigation of the mechanical unfolding of titin protein, for 
example, the restoring forces all show a sawtooth like pattern with a def- 
inite periodicity It reveals much information about the mechanism of the 
unfolding processes [331 HI] ■ The observed pattern, in addition to a fit of a 
worm-like chain model, has also been verified by the steered molecular dy- 
namics or Monte carlo simulations j3j3 HH] • Similar study has been extended 
to other systems 13711381133] . 

The same kind of experiments has also been performed by varying tem- 
perature. In the experiment carried out by Spider and Discher et al. |22] . 
the spectrin protein is chosen for the AFM study at different temperatures. 
Thousands of tip-to-surface contacts are conducted for a given temperature 
because of the statistical nature of the AFM measurement. The observed 
curve for the relation between the unfolding force and extension shows the 
similar sawtooth pattern for all temperatures. In addition, the tandem repeat 
unfolding events are more favored at lower temperature as demonstrated in 
the unfolding length histograms. Most striking is that the unfolding forces 
show a dramatically nonlinear decreasing relation as the temperature T ap- 
proaches the transition temperature T m . This is shown in Figure 3B of the 
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paper 

Similar behaviors regarding the force-temperature dependence have also 
been observed via either AFM or optical tweezers for the forced overstretch- 
ing transition for the individual double-stranded DNA molecules |2lH 1^1 EE 

EH- 

Some other interesting experiments which illustrate the effects from tem- 
perature on the microscopic structures of biomolecules have also been per- 
formed even though the detailed physical origins of the impacts (from either 
the electrons or the molecules) have not been specified jSHl EH EB EIH EH EE] ■ 
In a circular dichroism ( CD) spectra and high resolution NMR study, for in- 
stance, it shows that the secondary structure of the Alzheimer f3 (12-28) 
peptide is temperature-dependent with an extended left-handed 3i helix in- 
terconverting with a flexible random coil conformation [53]. Another example 
is related to the analysis of the temperature-dependent interaction of the pro- 
tein SshlOb with a DNA which influences the DNA topology [HU HUH EH! • The 
analysis from the heteronuclear NMR and site-directed mutagenesis indicates 
that the SshlOb exists as a dimer: T form and C form. Their ratio is deter- 
mined by the Leu 61 — Pro 62 peptide bond of the SshlOb which is sensitive to 
temperature. 

2.3 Temperature effects on intermolecular forces 

The analysis of the general issues related to the temperature effects on the 
microscopic structures has been most recently extended to the realm of in- 
termolecular forces. Since the usual intermolecular forces such as hydrogen 
bonds, van der Waals forces, ionic bonds, and hydrophobic interactions are 
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weak and typically of the order of 0.1 eV or 4-0 kT at the physiological 
temperature, the variation of temperature will thereby have a very strong 
influence over the strength of these forces. 

The first experimental investigation on the temperature-dependent inter- 
molecular forces is for the biotin-avidin system and by an AFM measure- 
ment jSlj. The biotin-avidin complex is a prototypical receptor and ligand 
system with the biotin binding strongly up to four avidin proteins jHHEOlinil 
IH21 EH] • They have an extremely high binding affinity, and therefore serves 
as a model system for various experimental examinations. In the experi- 
ment carried out in Beebe's group [HI], the receptor avidin is attached to the 
AFM tip and linked to the agarose bead functionalized with the biotin. The 
temperature of the entire AFM apparatus is varied at a range from 286 to 
"510K. In addition, the loading rate is kept very slow so that the thermal 
equilibrium for the biotin-avidin pairs is assumed. The forces expected to be 
determined is the rupture force Fi between the individual biotin-avidin pair 
which is defined as the maximum restoring force in an actual AFM ex- 
periment, however, the total adhesive force between the tip and substrate is 
measured. It is a sum of the finite number n of the interactions between each 
biotin and avidin pair. To extract the individual and average bond rupture, 
a statistical method has been developed in Beebe^s group [101 163] IB^ IBT] . 
They assume a Poisson distribution for the number n of the discrete rupture 
forces or linkages from multiple measurements, and have obtained the single 
force Fi at different temperatures. The result is shown in Figure 3 of the pa- 
per j64j . We see that the individual rupture force Fi for the biotin-avidin pair 
is decreased by about five-fold in strength when the temperature is increased 
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from 286 to 310K. 

To interpret the observed temperature impact on the biotin-avidin forces, 
Pee&e's group has performed a thermodynamic analysis Based on the 
simple models and arguments [SHUHH], they have come out an equation that 
connects the square of the single bond-rupture force Fi to the absolute tem- 
perature T as follows, 



where the k\, on d is the force constant of the individual biotin-avidin pair, and 
the time tr is the characteristic time needed to break n pairs of those forces. 
The is the energy required to remove the biotin from avidin's strongest 
binding site and the corresponding time is Tp. More details on this analysis 
can be found in the paper [64J. The relation between the square of the force 
Fi and the temperature T is also plotted as Figure 5 in that paper. Therefore, 
from the relation (1) and this figure, the information about the stiffness of 
the ligand and receptor bond and the critical binding energy, etc. can be 
obtained. Obviously, what we need at present is a microscopic theory which 
can account for all these relations and properties. 

2.4 Temperature effects on absorption spectra of elec- 
trons in solvents 

The structure and dynamics of the solute in solvents is one of the most 
important fields in chemistry since most of the chemical reactions occur in 
solution phases. In the meantime, it is also one of the most challenging 
areas in theoretical chemistry with many unsettled issues. The variation of 
temperature in the measurement of absorption spectra of solvated electron 
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in various solvents has proved to be a useful means for the understanding of 
the solvation processes [71-85]. 

There are several experimental techniques available for this type of in- 
quires with the pulse radiolysis being the one most commonly used. There 
are also several research groups conducting the similar experimental inves- 
tigations and obtaining the consistent results relating to the temperature 
effects on the optical absorption spectra of the solvated electron in solvents. 
In a recent experiment carried out in Katsumura's group, for example, the 
pulse radiolysis technique is employed to study the optical absorption spec- 
tra of the solvated electron in the ethylene glycol at different temperatures 
from 290 to 598-fT at a fixed pressure of 100 atm. In addition to the faster 
decay of absorptions, it is found that, their maximal positions shift to the red 
with the increasing temperature as shown in Figures 1 to 3 of the paper |S2] . 
This is in contrast to the situation for the electronic spectra of the polymers. 
They also point out the need to quantify the change of the density in the 
experiment in order to really understand the observed results. 

The same type of experiments has been extended to the examination of 
the optical absorption spectra for Ag° and Ag% in water by varying temper- 
ature, and similar results have been obtained [HTj . 




3 Theoretical Development 

Having presented four different types of experiments above, we can observe 
that the investigation of the temperature effects on the microscopic structures 
of molecules is a very interesting and sophisticated field. More need to be 
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probed and understood. Even though the experimental analysis has been for 
a long time, very limited number of the related theoretical works is available, 
especially at the first-principle level. In other words, the quantum chemistry 
at finite temperature is not a well-established field [13 EH El- 

It is true that the influence of temperature on the microscopic structures 
is a complicated phenomenon. There exists different functioning mecha- 
nisms. One consideration is that the variation of temperature, according to 
the Fermi-Dirac statistics, will change the thermal probability distribution 
of single-particle states for a free electron gas. It is expected that similar 
situation should occur for an interacting electron system, and therefore its 
microscopic structure will be correspondingly altered. Another considera- 
tion is that, for molecules or solids, the thermal excitation will cause the 
change of the time scales for the molecular motions. This will most likely 
bring about the transitions of the electronic states, and therefore lead to the 
breakdown of the Born-Oppenheimer approximation. The electron-phonon 
interaction is a fundamental topic in solid state physics and its temperature 
dependence is well-known. As a result, the variation of temperature will 
change the strength of the coupling between the electronic and molecular 
motions. Nevertheless, we tackle the issues pertaining to the temperature 
effects in a simpler way. We treat only an identical and interacting fermion 
system. Or we neglect the coupling of the electronic motion with those of the 
nucleus in the molecules or solids. We expect that some sort of the general 
conclusions will come out from this analysis. As a matter of fact, this is the 
approach commonly admitted in a non-adiabatic molecular dynamics study, 
in which purely solving the eigenequation for the electrons will provide the 
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reference states for the examination of the coupling between the electronic 
and nuclear motions of the molecules. 

In the following, we will present self-consistent eigenequations within the 
framework of the density operators in equilibrium statistical mechanics which 
decides the molecular orbitals at a given temperature. 

3.1 Hierarchy Bloch equations for reduced density op- 
erators in canonical ensemble 

We consider an identical and interacting iV-particle system. In a canonical 
ensemble, its iVth-order density operator takes the form 

D N = exp{-(3H N ), (2) 

and satisfies the Bloch equation [HU E2] 

- §pD» = H N D\ (3) 

where 

N N 

H N = J2h(i) + Y,9(h3), (4) 

i=l i<j 

is the Hamiltonian for the N particle system composed of one-particle op- 
erator h and two-body operator g. The f3 is the inverse of the product of 
Boltzmann constant kg and absolute temperature T. 

Since the Hamiltonian (4) can be written as a reduced two-body operator 
form, the second-order reduced density operator suffices to describe its N (> 
2) particle quantum states. A pth-order reduced density operator is generally 
defined by |H3 El 

TP = L P N (D N ), (5) 
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where L P N is the contraction operator acting on an iVth-order tensor in the 
iV-particle Hilbert space V . The trace of the D p gives the partition function, 



Tr(D p ) = Z(/3,V,N). 



(6) 



Rewrite the Hamiltonian in a form 

N p N N 

H N = H{+ £ fc(0+£ E 9(hj)+ E 9(i,j) 

j=p+l i=l i=p+l «<i(«>P+l) 



(7) 



where 

= EM*) + j), ( 8 ) 

i=l i<j 

and apply the contraction operator L P N on both sides of the Eq. (3), we 
develop an equation that the pth-order density operator satisfies [H] 



^D p = H p D p +(N-p)L p p+1 {h(p+l)D p+1 } + (N-p)L p p+1 



v_ 
Li=l 



J2g(hP + i)D p+1 



N -p 



^p+2 



g(p+l,p + 2)D p+2 



It provides a law according to which the reduced density operators vary in 
terms of the change of temperature. 



(9) 



3.2 Hierarchy Bloch equations for reduced density op- 
erators in grand canonical ensemble 

The above scheme for deducing the equations for the reduced operators can 
be readily extended to the case of a grand canonical ensemble. It is a more 
general one with a fluctuating particle number N. In this ensemble, the 
density operator is defined in the entire Fock space 



F = E ® v 

N=0 



N 
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and is written as the direct sum of the density operators D G (N) associated 
with the iV-particle Hilbert space V , 



where 



D G =J2 ®D G (N), 

N=0 



D G (N) = exp[-/3(H - //#)], 
= eM-PH), 



(10) 



(11) 



and 



H = H - fiN, (12) 

is called the grand Hamiltonian on V^. The form of the Hamiltonian H has 
been given by Eq. (4) and the \x is the chemical potential. The corresponding 
pth-order reduced density operator is therefore defined as 

' N ' 



U G 



N=p 



P 



L p n [Dg(N)} 



with the trace given by 



and 



Tr(D 



G) 



N 
P 



u Gi 



D 



o _ ■=>/ 



(13) 



(14) 



(15) 



. G -Z(P,H, V). 
The /x, V) is the grand partition function. 

In a similar manner, we can also derive the hierarchy equations that the 
reduced density operators in the grand canonical ensemble obey jH] 



.9 

dp 



HV D v + (p + 1)L p +i \h( p + i )jDP +ij + (p + 



.i=i 



p + 2 



L p p+2 g(p+l,p + 2)LP 



+2 



(16) 
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where 

m = Eh(i) + j:g(i,j), (17) 

i=l i<j 

and 

h(i) = h(i) - fi. (18) 

It gives us a law with which the reduced density operators in the grand 
canonical ensemble vary in temperature. 

3.3 Orbital approximation and Hartree-Fock equation 
at finite temperature 

The Eqs. (9) and (16) define a set of hierarchy equations that establish the 
relation among the reduced density operators D p , D p+1 , and D p+2 . They 
can be solved either in an exact scheme or by an approximate method. The 
previous investigation of iV electrons with an independent particle approx- 
imation to the Schrddinger equation for their pure states has lead to the 
Hartree — Fock equation for the molecular orbitals [95-101]. We thereby ex- 
pect that the same approximate scheme to the reduced Bloch equations (9) 
or (16), which hold for more general mixed states, will yield more generic 
eigenequations than the usual Hartree — Fock equation for the molecular or- 
bitals. 

We consider the case of a grand canonical ensemble. For p — 1, Eq. (16) 
reads 



-UpD 1 =HiD 1 +^^D 1 -l- D 1 hD 1 +2Ll [g(l, 2)D 2 } +3L\ [g(2,3)D 

(19) 

Under the orbital approximation, above second-order and third-order reduced 
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density operators for the electrons can be written as 



D A = D 1 A D 1 A D 1 /(D ^ 2 



(20) 



and 



D 2 = D 1 A D 1 / D° . (21) 

These are the special situations for the statement that a pth-order reduced 
density matrix can be expressed as a p-fold Grassmann product of its first- 
order reduced density matrices. With this approximation, the last two terms 
of Eq. (19) can be evaluated in a straightforward way as follows 



2L\ Ul,2)D 2 l ={J-K)D\ 



and 



where 



and 



3L\ 



<?(2,3)I7 



Tr{gD 2 ) 1 . 



~^D\J-K)D\ 



J = Tr<2 



D° D° 



g-D\2;2)] /D°, 



(22) 



(23) 



(24) 



K = Tro 



g-(2,3)-D 1 (2;2)\/D°, (25) 

are called the Coulomb and exchange operators, respectively. With (2, 3) 
being the exchange between the particle 2 and 3, the action of the K on the 
reduced density operator is 



K-D\Z;Z) = Tr 2 \g ■ (2,3) ■ D 1 (2; 2)1 / D Q ■ D l (3; 3) 



= Tr 9 



g-D\3;2)-D 1 (2;3)} / D°. 



(26) 



Substitution of Eqs. (22) and (23) into Eq. (19) results in the Block equation 
for the first-order reduced density matrix of the N interacting electrons under 
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the orbital approximation, 

d x„i (TrhD 1 Tr gD 2 \ . 1 . 

ni ^ e 1 ,,\ni i I i ^ I ni ni 



+ Z ^)B 1 -^C 1 (f'-^ 1 . (27) 



-_ D ' = (F-„) D ' +V D „ 
where 

F = h + J-K, (28) 

is called the FocA; operator at finite temperature. Redefine the normalized 
first-order reduced density operator 

p 1 = D 1 /D°, (29) 

we can simplify above equation into 

-■lip? = (30) 

Furthermore, from Eq. (30) and its conjugate, we get 

Fp 1 -p 1 F = 0, (31) 

which indicates that the Fock operator F and the first-order reduced density 
matrix p 1 commute. They are also Hermitian, and therefore have common 
eigenvectors These vectors are determined by the following eigen 

equation for the Fock operator, 

F\4>i) =ei|&). (32) 

It is the eigenequation for the molecular orbitals at finite temperature. 
The first-order reduced density operator is correspondingly expressed as 

P 1 = 5^(/3,/i,ei)|0i)(0i|, ( 33 ) 

i 
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where u(/3,[/,, e$) is the thermal probability that the orbital is found to be 
in the state {\4>i}} at finite temperature T. Substituting Eq. (33) into Eq. 
(30), we can obtain the equation this thermal probability u((3, fi, ei) satisfies, 

d 

- -gp^iP' ^> e i) = ( e * _ A^MA /i, Ci) - (ei - /i)^ 2 (/3, fl, ti). (34) 

Its solution has the same usual form of the Fermi — Dirac statistics for the 
free electron gas as follows, 

u(/3,fi,€i) = 1 + J,^) , (35) 
with the energy levels {e;} determined by Eq. (32). 



4 Discussion, Summary and Outlook 

In this paper, we have presented both experimental and theoretical develop- 
ments related to the temperature impacts on the microscopic structures and 
processes for the molecules. 

In the theoretical part, we have depicted the sets of hierarchy Block equa- 
tions for the reduced statistical density operators in both canonical and grand 
canonical ensembles for the identical fermion system with a two-body inter- 
action. We have solved the equations in the latter case under a single-orbital 
approximation and gained an energy eigenequation for the single-particle 
states. It is the extension of usual Hartree — Fock equation at absolute zero 
temperature to the one at any finite temperature. The average occupation 
number formula for each single-particle state is also obtained, which has the 
same analytical form as that for the free electron gas with the single-particle 
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state energy determined by the Hartree — Fock equation at finite tempera- 
ture (32). 

From Eqs. (24), (25) and (28), we see that the Coulomb operator J, 
the exchange operator K, and therefore the Fock operator F are both co- 
herent and incoherent superpositions of the single-particle states. They are 
all temperature-dependent through an incoherent superposition factor, the 
Fermi — Dime distribution, lj(/3, /i, 6j). Therefore, the mean force or the 
force field, and the corresponding microscopic structures are temperature- 
dependent. 

We have expounded the physics relating to the temperature effects on 
the electronic structures or processes for the molecules. These effects can be 
either at the single-electron level or of the iV-body excitation. This is very 
critical for our understanding and computation of the temperature influences 
over the molecular structures. From this analysis, for example, we can con- 
clude that the temperature should have a stronger effect on the transition 
states for the molecules. Accordingly, the change of the chemical reactivity 
might come from the alternation of the electronic states due to the variation 
of temperature. We expect that this sort of the changes involving configu- 
ration mixing will be very common for the systems we are discussing here. 
Therefore, it will be a very significant work to develop or test the correspond- 
ing multireference theory for the molecular orbitals or electronic structures 
at finite temperature [21 11U2| I1U3| 111)41 HU5j . More general or standard theo- 
retical explorations are expected [3]. 

On the experimental side, we have exposed four major fields of investiga- 
tion of chemical, biomolecular, or material importance, which demonstrate 
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the temperature impacts on their structures, spectra, or bond rupture forces. 

The complete determination of the geometric structures and electronic 
spectra of the polymeric molecules is a very difficult task. As has been 
stated in papers fOJ EI; there are many different elements contributing 
to the change of the spectra. At present, we focus on the examination of 
the effect from temperature. We have showed that it can alter both the 
shapes and positions of the absorption and other spectra for the PPV and 
its derivatives. As has been analyzed, the increase of temperature will bring 
about the excitation of the vibrational, rotational and liberal motions, which 
might also lead to the electronic transition. The Huang — Rhys parameter 
has been introduced to describe the strength of the coupling between the 
electronic ground- and excited-state geometries. Furthermore, it has been 
observed that this factor is an increasing function of temperature [HI El ED] . 
Obviously, a more detailed analysis of the electronic structure, excitation 
and spectroscopic signature at the first-principle level, which includes the 
temperature-dependent force field, is anticipated. 

Temperature has proved to be a big player in both experimental and 
theoretical study of the structure and dynamics of biomolecules including 
their folding or unfolding. At a first glance, the energy gap between the 
HOMO and LUMO for the biomolecules should be small or comparable to 
the Boltzmann thermal energy ksT because of their very large molecular size. 
Therefore a change of temperature should have a strong influence over their 
electronic states, and consequently, the energy landscape and the related dy- 
namics including the folding and unfolding. The experimental investigation 
with the AFM and other techniques of the temperature effect on the shift of 
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their unfolding pathways might have verified this type of thermal deforma- 
tion of the potential energy landscape fK)\ l2*Tj 12*2*]. This is in contrast to the 
tilt and deformation of the energy landscape including its transition states 
for the biomolecules resulting from the applied mechanical forces |7Uj . 

Unfolding proteins by temperature is not just one of the classical experi- 
mental techniques for the analysis of the structure, dynamics and energetics 
of the biomolecules. It has also been utilized, for example, in the molec- 
ular dynamic simulation to study the structure of the transition states of 
CI2 in water at two different temperatures: 298 K and 498 K [23] ■ The 
later high temperature is required in order to destabilize the native state for 
monitoring the unfolding as done in the real experiments. In another re- 
cent molecular dynamics simulation [23], Karplus's group has compared the 
temperature-induced unfolding with the force-stretching unfolding for two (3- 
sandwich proteins and two a-helical proteins. They have found that there are 
significant differences in the unfolding pathways from two approaches. Nev- 
ertheless, in order to get more reliable results, the temperature-dependent 
force fields need to be developed and included in the molecular dynamics 
simulations. This is also the case in the theoretical investigation of protein 
folding since an accurate simulation of protein folding pathways requires bet- 
ter stochastic or temperature-dependent potentials which have become the 
bottleneck in structure prediction [23 I2H1 123 EE [2Hj • From structural points 
of view, the variation of temperature leads to the change of the mean force 
or the energy landscape, and therefore provides a vast variety of possibilities, 
for instance, in the protein design and engineering. 

The intermolecular forces are ubiquitous in nature. They are extremely 
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important for the biological systems and for the existence of life. The inter- 
molecular forces have the specificity which is responsible for the molecular 
recognition between the receptor and ligand, the antibody and antigen, and 
complementary strands of DNA, and therefore for the regulation of complex 
organization of life |1(J6| . For these reasons, the experiment carried out in 
Beebe's group has an immediate significance. It has demonstrated that tem- 
perature can be an important factor for changing the specificity of the inter- 
molecular forces and therefore the function of life [H3j . Nevertheless, how the 
charge redistribution occurs due to the variation of temperature has not been 
interpreted, and a microscopic theory for quantifying the temperature influ- 
ence over the intermolecular forces is still lacking. Since the delicate study 
of the intermolecular forces provides the insight into complex mechanisms of 
ligand-receptor binding and unbinding processes or pathways, a paramount 
future research is to establish the links between the intermolecular forces and 
the temperature within quantum many-body theory. 

The theoretical exploration of the temperature effects on the optical ab- 
sorption spectra of solvated electrons is still in very early stage and few 
published works are available [BI3 EH EH] • One of the earliest inquires 
by Jortner used a cavity model to simulate the solvated electron where the 
electron is confined to the cavity surrounded by the dielectric continuum sol- 
vent [SSI- However, his investigation is not of fully microscopic in nature 
since he assumed a temperature dependence of phenomenological dielectric 
constants which were obtained from the available experimental data. In ad- 
dition, the model used is too simplified and, for instance, it neglects the 
intrinsic structure of solvent molecules. There are a few recent examinations 
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on the temperature influences over the absorption spectra of the solvated 
electrons. They all cannot catch the full features of the experimental obser- 
vations. One reason is that the physical nature for the process is not totally 
understood which might leads to incorrect models employed for the simu- 
lation. The other is to utilize the crude models which might have omitted 
some important physical effects. For example, in an analysis by Brodsky and 
Tsarevsky [S3], they have concluded a temperature-dependence relation for 
the spectra which is, however, in contradiction with the experimental find- 
ings at high temperature. The quantum path-integral molecular dynamics 
simulation cannot produce those temperature-dependence relations observed 
in the experiments [S3]. In a recent quantum-classical molecular-dynamics 
study by Nicolas et al, even though the temperature-dependent features of 
optical absorption spectra for the solvated electron in water have been recov- 
ered jHEl, however, they claim that the red shifts of absorption spectra with 
the increasing temperature observed in both experiments and calculations 
are due to the density effect instead of temperature. This might cast the 
doubt of usefulness of our present theoretical work in this area. However, 
after examining their work, we observe that they actually have not included 
any temperature effects on the electron in their theoretical model. These 
effects might be either from the Fermi — Dirac distribution for individual 
electrons or due to the electronic excitation caused by the thermal excitation 
of the solvent, as we have discussed previously and in paper jUJ. Obviously, 
much finer theoretical works or more experimental investigations in this area 
are expected to resolve this dispute. 

In addition to the systems discussed above, there are many other types 
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which show the temperature impacts on their microscopic structures. Either 
theoretical or experimental works have been done or are in progress. Exam- 
ples include the study of the temperature dependence of the Coulumb gap and 
the density of states for the Coulomb glass, the experimental investigation 
of the temperature effects on the band-edge transition of ZnCdBeSe, and 
the theoretical description of the influence from temperature on the polaron 
band narrowing in the oligo-acene crystals |1U7| 1108] 109 . 

To sum up, the quantum chemistry at finite temperature is a new and 
exciting field. With the combination of the techniques from modern quantum 
chemistry with those developed in statistical or solid state physics, it will 
provide us with a myriad number of opportunities for the exploration [3J 
ITMITUSllTTn] 1 
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